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Abstract 

Genetically identical cell populations exhibit considerable intercellular variation in the level of a given protein or mRNA. 
Both intrinsic and extrinsic sources of noise drive this variability in gene expression. More specifically, extrinsic noise is the 
expression variability that arises from cell-to-cell differences in cell-specific factors such as enzyme levels, cell size and cell 
cycle stage. In contrast, intrinsic noise is the expression variability that is not accounted for by extrinsic noise, and typically 
arises from the inherent stochastic nature of biochemical processes. Two-color reporter experiments are employed to 
decompose expression variability into its intrinsic and extrinsic noise components. Analytical formulas for intrinsic and 
extrinsic noise are derived for a class of stochastic gene expression models, where variations in cell-specific factors cause 
fluctuations in model parameters, in particular, transcription and/or translation rate fluctuations. Assuming mRNA 
production occurs in random bursts, transcription rate is represented by either the burst frequency (how often the bursts 
occur) or the burst size (number of mRNAs produced in each burst). Our analysis shows that fluctuations in the transcription 
burst frequency enhance extrinsic noise but do not affect the intrinsic noise. On the contrary, fluctuations in the 
transcription burst size or mRNA translation rate dramatically increase both intrinsic and extrinsic noise components. 
Interestingly, simultaneous fluctuations in transcription and translation rates arising from randomness in ATP abundance 
can decrease intrinsic noise measured in a two-color reporter assay. Finally, we discuss how these formulas can be 
combined with single-cell gene expression data from two-color reporter experiments for estimating model parameters. 
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Introduction 

Genetically identical cell populations exposed to same extracel- 
lular environment exhibit considerable variability in gene expres- 
sion [1-5]. This variation in the level of a given protein is often 
referred to as gene expression noise. Increasing evidence suggests that 
noise plays important functional roles in many cellular processes. 
For example, tight control of expression noise is vital for optimal 
functioning of housekeeping proteins [6-8] , and diverse diseased 
states have been attributed to an elevated expression noise [9-1 1]. 
Not surprisingly, genes actively use different regulatory mechanism 
to reduce stochastic fluctuations in protein levels [12-22,22-25]. 
Expression noise is also exploited to drive genetically identical cells 
to different cell-fates [26-31], and to buffer cellular populations 
from hostile changes in the environment [27,32—34]. 

Gene expression noise can be decomposed into intrinsic and 
extrinsic noise [35-37]. More specifically, intrinsic noise is the 
protein variability that arises from the inherent stochastic nature of 
biochemical reactions associated with transcription, translation, 
mRNA and protein degradation. Given that many mRNA species 
are present at low copy numbers inside cells, random birth and 
death of individual mRNA transcripts generates considerable 
intrinsic noise [38-41]. Let ^ be any cell-specific factor (such as 
cell cycle stage, abundance of RNA polymerases/ribosomes, 



cellular environment, etc.) that affects expression of a given gene. 
Then, cell-to-cell differences in Z wu l create intercellular 
variability in gene expression, that is referred to as extrinsic noise 
Variations in £ induce fluctuations in model parameters (such as 
the transcription and translation rate), and extrinsic noise can be 
effectively quantified through analysis of deterministic gene 
expression models with corresponding parameter fluctuations [42] . 

We define intrinsic and extrinsic noise in the context of a two- 
color experiment, where the gene of interest is duplicated inside 
the cell (Figure 1). Consider two identical copies of a promoter that 
express two different reporter proteins Pi and Pi- Let P\(f) and 
P2(t) denote the level of these proteins at time t inside the cell. 
Since cell-specific factor Z is common to both copies of the gene, 
cell-to-cell variations in Z will make P\(t) and P2(t) correlated. 
The contribution of Z to expression noise is quantified via the 
extrinsic noise defined as 

Extrinsic Noise = <^>-<Pi><^ 

<pO <>2> 

and is related to the covariance between reporter levels. If reporter 
levels are perfectly correlated, and assuming (p\} = (pi), 

(P\~) = <P 2 2>, 
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Extrinsic Noise = Total Noise : 



<p]}-<po 2 

<pi> 2 



(2) 



which is the total noise in protein level measured by its coefficient 
of variation squared. Intrinsic noise is the protein variability that is 
not accounted for by extrinsic noise, and is defined as 



Intrinsic Noise = Total Noise — Extrinsic Noise 

= <I> 2 > - <I>iP2> 
<Pi> 2 



(3) 



In summary, a two-color assay can be used to decompose the 
total protein noise level into intrinsic and extrinsic noise 
components, computed via (1) and (3), respectively. 

Analytical formulas for intrinsic and extrinsic noise are derived 
for a class of stochastic gene expression models with fluctuations in 
the transcription or translation rate. Assuming mRNA production 
occurs in random bursts, transcription rate is represented by either 
the burst frequency (how often the bursts occur) or the burst size 
(number of mRNAs produced in each burst). Our results show that 
fluctuations in the transcription burst frequency enhance extrinsic 
noise but do not affect the intrinsic expression noise. However, 
fluctuations in the transcriptional burst size or mRNA translation 
rate increase both intrinsic and extrinsic noise. A recent study has 
implicated fluctuations in ATP levels as a major driver of gene 
expression variability [43]. Since ATP affects both transcription 
and translation, simultaneous fluctuations in multiple model 
parameters is investigated. Interestingly, simultaneous fluctuations 
in the transcription and translation rates decrease intrinsic noise in 
certain parameter regimes. Finally, usefulness of these formulas in 
interpreting two-color reporter experiments and estimating model 
parameters is discussed. 

Gene Expression with Constant Parameters 

We begin by introducing the standard stochastic gene expres- 
sion model [44—47], where all model parameters are fixed, and 



expression variability arises due to the stochastic nature of 
transcription and translation processes. 

Model Formulation 

Transcription has been shown to occurs in "bursts" with each 
burst producing multiple mRNA copies [48-53]. Assume mRNAs 
are produced in bursts of size B m that occur at a rate k m . We refer 
to k m and B m as the transcriptional burst frequency and burst size, 
respectively. Consistent with measurements [50], B m is assumed to 
be a geometrically distributed random variable with probability 
distribution 

Probability{S m = ;} = a, t = (1 -sfs, 0<s<\, ; = {0,1,2, . . .} (4) 

and mean burst size <-B»,> : =(\—s)/s. Proteins are produced 
from each mRNA at a translation rate k p . Finally, mRNAs and 
proteins degrade at constant rates y m and y p , respectively. The 
stochastic model considers transcription, translation and degrada- 
tion as probabilistic events that occur at exponentially-distributed 
time intervals [54,55]. Moreover, whenever a particular event 
occurs, the mRNA and protein population count is reset 
accordingly. Let m(t) and p(t) denote the number of molecules 
of the mRNA and protein at time t, respectively. Then, the reset in 
m(t) and p(t) for different events is shown in the second column of 
the table in Figure 2. The third column lists the propensity 
functions f(m,p) which determine how often an event occurs. In 
particular, the probability that a particular event will occur in the 
next infinitesimal time interval (t,t + dt\ is given by f(m,p)dt. 

Computation of Intrinsic Noise 

It is relatively straight forward to derive differential equations 
describing the time evolution of the different statistical moments of 
the mRNA and protein count. For the above model, the time- 
derivative of the expected value of any differentiable function 
<p(in,p) is given by 

^^>=<^^( W1 ,rtx/(^)>, (5) 

Events 

where Acp(m,p) is the change in cp when an event occurs, f(m,p) is 



Protein (~) A Protein 

t t 




Extrinsic Factor Fluctuations 



Total Noise = Variance in protein levels 
Mean protein level 2 

Extrinsic Noise = Covariance in protein levels intrinsic Noise = Total Noise - Extrinsic Noise 
Mean protein level 2 



Figure 1. Decomposing gene expression variability into extrinsic and intrinsic noise using a two-color reporter assay. Two identical 
copies of a promoter express two different reporter proteins. Correlation in reporter levels is a measure of extrinsic noise that arises from cell-to-cell 
differences in shared cellular factors. Intrinsic noise is the protein variability that is not accounted for by extrinsic noise, and typically originates from 
the inherent stochastic nature of biochemical processes. 
doi:10.1371/journal.pone.0084301.g001 
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Figure 2. Model formulation. Schematic of the gene expression model (left). The stochastic model consists of four events that occur randomly at 
exponentially-distributed time intervals. Discrete changes in the mRNA (m(t)) and protein [p(t)) population count for different events are shown in 
the second column of the table. Third column lists the event propensity function that determines how often an event occurs. 
doi:1 0.1 371 /journal, pone.0084301 .g002 



the event propensity function, and <.) represents the expected 
value [56,57]. Using the resets and propensity functions in Figure 2 
this corresponds to 



used as a metric for quantifying the extent of variability/ noise in 
protein copy numbers. From the steady-state protein variance and 
mean we obtain 



d<<p(m,p)} _ 
dt 

<y m m W{m -Up)- (p{m,p)] + k p m[<p(m,p + 1) - <p{m,p)] > 
+ <.VpPl9(mj> - 1) - <p(m,p)] > 

00 

+ <^ k m ai[g>(m + i,p) - <p(m,p)}}. 



(6) 



CVi 



fixed ' 



(Bjy + (B,„y y„ | J_ 
2<5„,><m> y p +y m <j>y 



(9) 



which represents the total intrinsic noise in protein level for fixed 
parameters. As B m is geometrically distributed, (B 2 ,y = 
2(B m y 2 + (B m }, and (9) reduces to 



Choosing (p(m,p) as m,p,m 2 ,p 2 and mp in the above equation 
yields 



=k m iB m }-y„<jn}, =k p (m) -y p <p} (7a) 



=k l „(B 2 „y + y, n <,m> + 2<.B m )k m (my-2y m (m 2 y (7b) 



d<p 2 y 



dt 



-- k P <m> + y p <p y + 2k p (mp)- 2y p (j? 2 > ( 7c) 



d(mpy 
dt 



--k p <m 2 }+ (B m yk m <j>y-y p (mpy-y m {mpy. (7d) 



Setting the left-hand-side of (7) to zero and solving for the 
moments results in the following steady-state mean protein and 
mRNA levels 



— — k m (B m y k p (m) 



(8) 



where (,B m y is the mean transcriptional burst size and <(.) 
represents the steady-state expected value. As done in previous 
studies of intrinsic and extrinsic noise [35,36,58], the steady-state 
coefficient of variation squared (variance divided by mean squared) is 



cv- 



2 _<£„,> + 1 y p 

fixed 



(my y p +y„ 



1 

w 



(10) 



The first term on the right-hand-side of ( 1 0) represents the noise in 
mRNA copy numbers that is transmitted to the protein level 
[15,46]. The second term is the Poissonian noise arising from 
random birth-death of protein molecules. Next, the noise 
additional to (10) that comes from fluctuations in individual model 
parameters (such as k m , (B ln y and k p ) is quantified. 

Transcription Burst Frequency Fluctuations 

Consider a cell-specific factor Z at the transcriptional level (such 
as a transcription factor). Then, fluctuations in Z can either affect 
the transcriptional frequency k m or burst size B m in the model. 
The former case of burst frequency fluctuations is considered first. 

Modeling Parameter Fluctuations 

Let z(t) denote the level of a cellular factor Z inside the cell at 
time t. Fluctuations in z(t) are modeled through a simple birth- 
death process with probabilities of formation and degradation in 
the infinitesimal time interval (1,1 + dt] given by 

Probability {z(t + dt) = z(t)+\} = k z dt (11a) 

Probability {z(t + dt) = z( t) - 1 } = y z z(t)dt, (lib) 

where k z and y z represent the production and degradation rate of 
Z, respectively. For the process described in (1 1), the steady-state 
mean, coefficient of variation squared CV 2 and the auto- 
correlation function .R-(t) are given by 
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<z>=^, 0-f=JL 7? z (T)=exp(-y z T). (12) 

Yz <z> 

Thus by changing k z and y z , both the extent and time-scale of 
fluctuations in z(t) can be independently modulated. Note the 
inverse relationship between <z> and CV 2 implies Poisson statistics. 
Fluctuations in Z are incorporated in the model by assuming that 
the transcription burst frequency is no longer a constant but given 
by k m z(t)/ \z~), making it a random process with mean k m and 
coefficient of variation squared CV 2 . Throughout this manuscript, 
CV 2 represents the extent of parameter fluctuations. Since Z 
similarly affects expression of both copies of the gene in a two-color 
assay, fluctuations in z(i) make reporter levels correlated in Figure 1 
and induce extrinsic noise. 

Computation of Total Noise 

The stochastic model consists of six birth-death events that 
change cellular factor, mRNA and protein copy numbers by 
integer amounts. Using the propensity functions in Figure 2 and 
(11) in (5) we obtain 



d((p(m,p,z)) 
dt 



= (y z z[<p(m,p,z- \)-q>(m,p,z)\y 

+ (k z [q>{m,p,z+ Y)-q>(m,p,z)]y 
+ <.Y P P[<P(m,p - 1 ,z) - <p{m,p,z)] > 
+ (k p m[(p(m,p + l,z) — <p(m,p,z)] > 
+ <y m m[q>{m - 1 ,p,z) - <p(m,p,z)] > 

k z 00 
w i=0 



(13) 



^^ = ^<m 2 >+/c„,<5„,><pz>/<z>-y p <m/ ) >-y m <^>(14e) 



<3?<OTZ> 

dt 



--kAm} + k m <:B m y(,z 2 y/(zy-y m (,mzy-y z <:mzy (14f) 



dipzy 



=k z <j>y+k p <wz> - y p <j>zy-y z (pzy, ( 14g) 



which yield the steady-state variability in protein level as 



CV 2 

burst— freq. 



<B m y+i 



■ + : 



1 



+ CV; 



( m y y p +y m <j>y 

2 YmVpiVm+Yp+Yz) 
(Ym+Yp)(Ym+Y z )(Y z + YpY 



(15) 



The first two terms on the right-hand-side of (15) represent the 
noise level with fixed parameters (Eq. (10)). The third term is the 
additional noise due to burst frequency fluctuations. Next, (15) is 
decomposed into intrinsic and extrinsic noise components as 
measured by the two-color reporter assay (Figure 1). 

Computation of Intrinsic and Extrinsic Noise 

Extrinsic noise can be approximated by the coefficient of 
variation squared of the protein level in a deterministic gene 
expression model with corresponding parameter fluctuations [42] . 
The deterministic counterpart to the stochastic model is the set of 
ordinary differential equations 



dm(i) k m (B m yz{i) 



dt 



<z> 



-Y m m(t) 



(16a) 



for any diflerentiable function <p(m,p,z). Appropriate choices of 
(p(m,p,z) result in 



dp(t) 



-k p m{t)-y p p{t) 



(16b) 



d ^ =k z -y z (zy, d ^J^ =k m (B m yz/ (zy-y m (my, 
=k p (my-y p (py 



(14a) 



d<z 2 y 
dt 

d(m 2 y 

~^r '~ ^" ' (14c) 

y m (my + 2k m (B m y(mzy/{zy-2y m (m 2 y 



-- k z + y z <z> + 2k z <z> - 2y,<z 2 > (14b) 



-k m (B 2 n y(zy/(zy + 



driven by the stochastic process z(t) defined in (1 1). For this hybrid 
model, where some states are continuous and other are discrete, 
the time derivate of {(p(m,p,z)y is given by (see Theorem 1 in [59]) 



dt 



-(y z z[(p(m,p,z- \)-<p{m,p,z)]y 

+ (k z [(p(m,p,z+l)-<p(m,p,z)]} (17) 

B(p(m,p,z) fk m zB m \ df{m,p,z) , . 
+ < ^7. ( -^T- - Ym m + T~„ ( k P m - YpP) > 



dm V <z> 

and leads to moment dynamics identical to (14) except for 



d(m 2 y 
dt 



= 2k m (B ln y(mzyi(zy-2y m (m 2 y (18a) 



d<p 2 y 

dt 



= k p (my + y p <py + 2k p On P y-2y p (p 2 y ( 1 4d ) 



d ^=2k p (m P y-2y p (p 2 y. 



(17) 
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Quantification of protein noise level from (14) (with (14c)-(14d) 
replaced by (18a)-(18b)) gives the extrinsic noise, which is 
subtracted from (15) for the intrinsic noise. This analysis results in 



where a, is given by (21). Equation (24) yields moment dynamics 
identical to (14) except for the time derivative of (rn 2 (t)y. For 
(p(m,p,z) = m 2 , 



Total noise = CV 2 hurst _ freq 

= Intrinsic noise + Extrinsic noise 



(19a) 



d(m 2 ) 
dt 



: 7m (m>- 2y m (m 2 >+ k m <^ aj 2 >+ 2k m (m ^ a,->(23) 



1 = 0 



i = 0 



Intrinsic noise = CV, 



fixed ' 



<B,„y+i Yp | i 

<m> Vp + Vm <P> 



, • ■ y--iT^2 YmYp^Ym Yp Yz) 

Extrinsic noise = CV, '— '— . 

(Ym + YpKYm + YzAYz + Yp) 



(19b) 



(19c) 



Using the fact that for a geometric distribution 



a,; 2 =2i y^q,-; ) + a,- 



1=0 \i=0 

and (21), (23) is written as 



(24) 



As expected, extrinsic noise increases with extent of parameter 
fluctuations CV 2 . On the contrary, intrinsic noise is independent 
of CV 2 and is equal to CV 2 ixed . An important limit considered 
previously is the case where parameter values (in this case 
transcription burst frequency) are drawn from a static distribution 
[36]. In our model, this corresponds to a scenario where the time- 
scale of fluctuations in z(t) are slow compared to mRNA/protein 
turnover rates. When Yz^Ym'Ypt Eq. (19c) reduces to 
'Extrinsic noise = CV 2 , and this result is consistent with previous 
calculations of extrinsic noise for parameter values drawn from a 
static distribution (see Eq. 25 in [36]). 

Transcription Burst Size Fluctuations 

Consider an alternative scenario of a fixed transcription burst 
frequency but varying burst size. Assume mRNAs are produced in 
geometrically distributed bursts with mean <_B„,>z(?)/<z>, where z(t) 
is the level of the cellular factor inside the cell at time t. This implies 

Probability{£ m = i} = o: i = (l-.?(0)'.<0, i = {0,l,2, . . .}, (20) 
and mean burst size 



\-s(t) _(B m }z(t) 
<t) ~ <z> 



>s(t) = 



1 



<B m )z(t) ■ 



1 + 



(21) 



dt = Ym <m>~ 2y m <m > + 



(2CB m ><z 2 > + <z><z> + 2<mz><z» . 



(25) 



Steady-state analysis of (14) (with (14c) replaced by (25)) results in 



cvl 

^ burst— size 



(B m Xl + CV 2 )+l Y P 

On} Yp + Yn 

1 , CV 2 YmY P (y m +Y P +Yz) 

<F> Z (Ym + Yp)(Ym + Yz)(Yz + YpY 



■ + 



(26) 



the total protein noise level for transcriptional burst size 
fluctuations. As expected when CV 2 = 0 (no parameter fluctua- 
tions) (26) reduces to (10). Comparison of (26) with (15) reveals that 
for a given CV 2 , burst size fluctuations generates larger variability 
in protein level than burst frequency fluctuations. 

Computation of Intrinsic and Extrinsic Noise 

For burst size fluctuations, the deterministic model used for 
quantifying extrinsic noise will be identical to (16). Since both 
transcriptional burst size and frequency appear together, replacing 
£ m z(r)/<z> with k m , and (B m } with (B m }z(t)/(z} in (16) does 
not alter the model. Thus, extrinsic noise is same irrespective of 
whether fluctuations are in the transcriptional burst size or 
frequency. Using (19c) and (26) 



Computation of Total Noise 

Time derivative of statistical moments is obtained from 



Total noise=CK / 2 „ ra( _ SIZ( ,=Intrinsic noise+Extrinsic noise (27a) 



dt 



+ {k : [q>(m,p,z + 1 ) - (p(m,p,z)\ > 

+ <Y P P{<P(m,p - 1 ,z) - <p(m,p,z)] + k p m[tp(m,p + l,z) - (p{m,p,z)\ >p 3 ) 
+ <Ym m l<P( m ~ 1 >P> Z ) ~ ?( m ,P^)\ > 

CO 

+ <C^2 k m a.j[(p{m + i,p,z) - <p{m,p,z)} >, 



Intrinsic noise= 



{my Y p + Ym <py fixedy 



Extrinsic noise = CV: 



2 YmYp (Ym + Yp + Yz) 



(Ym+ YpKYm +Yz)(Yz + Yp) 



(27c) 



In contrast to (19), intrinsic noise linearly increases with CV 2 for 
burst size fluctuations (Figure 3). 
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Figure 3. Gene expression variability for individual-parameter fluctuations. Intrinsic and extrinsic noise measured in two-color assay as a 
function of CV} (extent of parameter fluctuations) for fluctuations in the transcription burst frequency (left), transcription burst size (middle) and 
mRNA translation rate (right). Intrinsic noise is independent of CV} for transcription burst frequency fluctuations. However, for transcription burst 
size or translation rate fluctuations, intrinsic noise increases with CK 2 . Extrinsic noise always increases with CK_ 2 and is the largest for translation rate 
fluctuations. 

doi:1 0.1 371 /journal.pone.0084301 .g003 



Translation Rate Fluctuations 

Next, we consider mRNA translation rate fluctuations and set it 
equal to k p z(t)/ '(z). From Figure 2, this implies that the 
propensity function for the translational event is now nonlinear 
and given by kpZ(t)m(f) j \z} . Since mRNA production is no 
longer dependent on Z, z(t) and m(t) are independent random 
processes. 

Computation of Total Noise 

Statistical moments of z(t),m(t),p(t) are obtained from (13) with 
A: m z(f)/<z> replaced by k m , and k p replaced by k p z(t)/(zy. Using 
the fact that z(t) and m(t) are independent yields 



k--y-(z}, ——- =k m {B m y-y m (m), 



dl 
d<p> 
dt 



dt 



--k p (mXz>/<z>-V P <P> 



(28a) 



d(z 2 ) 
dt 



--k z +y z (z} + 2k-_ <z> - 2y z <z 2 > (28b) 



^ 2 >=^<^ m >< Z >+y m <m> + 
2k m (B m }(mz > - 2y m <m 2 > 



d<p 2 } 



-k p <,mXz>/<z>+Y p <P> 



+ 2k p (rnpz}/ <z> - 2y p <j> 2 } 



d -^=k p<m iy ( zym + 

k m (B,„ y<j?zy/(z}-y p (mp} - y m (mp) 



(28c) 



(28d) 



(28e) 
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d(pz) 
dt 



= k z < P y + k I ,(mXz> 2 /(z)-y p <pz>-y z (pz). (28f) 



Total noise = CVl anslation _ mte 

= Intrinsic noise + Extrinsic noise 



(34a) 



Note that the moment dynamics is not closed, in the sense that, 
the time derivative of the second order moments (p 2 {t)y depends 
on the third order moment (rn(i)p(t)z(t)y . This phenomenon 
occurs due to nonlinear propensity functions and typically closure 
methods are needed to solve for the moments [56,57]. The 
independence of z(f) and m(t) is exploited for moment closure. 
More specifically, 



Intrinsic noise = 



<B m y+\ 

<OT> 



1 + 



Yv 



cv?(y m +v p ) 

Ym+Yp+V Z I Ym+Yp ' <F> 



+ : 



(34b) 



dijnpz) 
dt 



--k z (mp}+k m (B m }<j>zy + 



(29) 



k p (m 2 z 2 y/ (zy-y z (mpzy-y p (mpzy-y m (mpzy 
which is dependent on the fourth order moment <wi 2 z 2 >. As 

<m 2 z 2 > = <m 2 ><z 2 >, (30) 

equations (28)-(30) form a closed system of equations that yield 
total variability in protein level as 



CV; 



translation — rate 



+ 



<B m y + i 



Yp + Y z 
. CV}y p 



(31) 



■ + ■ 



{niy \Y m + Y P Y m + Y P + Yz 



+ ■ 



Computation of Intrinsic and Extrinsic Noise 

Strategy for decomposing (31) into its intrinsic/extrinsic 
components is similar to previous sections: extrinsic noise is first 
computed from a deterministic model and then subtracted from 
(31) for the intrinsic noise. Consider the differential equation 
model 



dm(t) 
dt 



--k„(B m y-y m m(t) 



=k p m{t)z(t)/(zy-y p p(t). 



(32a) 



(32b) 



with translation rate fluctuations. Replacing k m z(t)/(zy by k m , 
and k p by k p z{t)/(zy in (17), we obtain moment dynamics 
identical to (28) except for 



d(m 2 y 
dt 



d<p 2 y 
dt 



= 2k m (B m y(mzy-2y m (m 2 y (33a) 



= 2k p (mpzyi(zy-2y p (p 2 y. 



(33b) 



Steady-state analysis of (28)-(30) (with (28c)-(28d) replaced by 
(33a)-(33b)) yields 



Extrinsic noise = 



_ CV 2 y p 
Y P +Y Z ' 



(34c) 



As in (27), fluctuations in the translation rate enhance both 
intrinsic and extrinsic noise (Figure 3). 

Simultaneous Model Parameter Fluctuations 

Previous sections focused on expression variability generated by 
fluctuations in individual parameters. However, stochasticity in the 
abundance of certain cellular factors (such as ATP) can 
simultaneously affect both transcription and translation. Motivated 
by this scenario, we investigate how perfectly correlated fluctua- 
tions in the transcription rate (measured by either the transcrip- 
tional burst frequency or burst size) and translation rate affect 
intrinsic and extrinsic noise. 

Transcription Burst Frequency and Translation Rate 
Fluctuations 

Assume transcriptional bursts occur at a rate k m z(t)j <(z)> with a 
geometrically distributed burst size independent of z(t) and given 
by (4). Each mRNA produces proteins at a rate k p z(i)j(z.y, which 
is perfectly correlated with burst frequency. Let 

LL=[(zy,(my,<jjy,(z 2 yxm 2 y,<j, 2 y,(mzy,<jjzy,(m P y] T (35) 

be a vector containing all the first and second order moments of 
the population counts. Then, using (13) with k p replaced by 
k p z(t)/(zy, time evolution of p can be compactly represented as 

d -^ = a x +A x H+Blfi, Ji = [<mpzy,<j>z 2 >,<m 2 z>,<mz 2 >,<z 3 >] r (36) 

where vector d\, matrices A\, B\ depend on model parameters 
and p is a vector of third order moments. As one would expect, 
nonlinear propensity function for the translation event leads to 
unclosed moment dynamics. It turns out that incorporating certain 
higher order moments in p can close moment equations. More 
specifically, the time derivative of 



£=[/,^,<™ 2 z 2 >,<mz 3 >,<z 4 >] r 
is closed and is given by 

dp. 



dt 



= 03+^3/' 



(37) 



(38) 



for some vector 03 and matrix At,. Steady-state analysis of (38) 
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O 

o o 
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0 0.5 1 

Extent of Parameter Fluctuations 



Figure 4. Gene expression variability for multiple-parameter fluctuations. Intrinsic noise measured in two-color assay as a function of CV: 
(extent of parameter fluctuations) for simultaneous fluctuations in the transcription burst frequency/translation rate (left), and transcription burst 
size/translation rate (right). The latter case generates larger intrinsic noise and also yields different qualitative trends compared to burst frequency/ 
translation rate fluctuations. Depending on parameter regimes, intrinsic noise can increase, decreases or change non-monotonically with CK Z 2 . High, 
medium, low protein populations correspond to an average of 300, 30 and 1 0 protein copies per cell, respectively. Other model parameters taken as 
mRNA half-life = 2 hours, protein half-life = time-scale of parameter fluctuations = 10 hours, mean transcriptional burst size = 10 and mean mRNA copy 
number per cell = 50. 
doi:1 0.1 371 /journal.pone.0084301 .g004 



results in an exact analytical formula for the total steady-state 
protein noise level. In previous sections (individual parameter 
fluctuations), average protein copy number was invariant of CK? 
and given by (8). However, simultaneous transcription/translation 
rate fluctuations enhance mean protein level from (8) to 



</>> = 



Ynu'n 



CVh 



-7z 



(39) 



To resolve total noise into its intrinsic /extrinsic components the 
following deterministic model is used 



be listed here but are illustrated in Figure 4. Interestingly, 
simultaneous fluctuations in the burst frequency and translation 
rate can either increase or decrease intrinsic noise depending on 
model parameters. 

To further elucidate the relationship between intrinsic noise and 
CVl, the case of slow fluctuations in z(t) compared to mRNA/ 
protein turnover rates (i.e., y z «y m ,y p ) is considered. In this case 
noise expressions reduce to 



Intrinsic noise = 



_<£,„)+ 1 



CVt 



(l + CV^YJY, 



+=L(41a) 



dm(t) k m (B„,yz(t) 
dt ~ <?> 



dp(t) _ k p m(t)z(t) 



dt 



<z> 



<p- 



,p(t)- 



(40a) 



(40b) 



For (40), the moment generator equation is obtained by 
replacing k p with k p z(i)j '<z> in (17). Performing an identical 
analysis as (36)-(38) for the hybrid model (40) yields the extrinsic 
noise, which is subtracted from the total noise to obtain the 
intrinsic noise. Unfortunately, these expressions are too complex to 



CK?(4 + 6CF. 2 + CK. 4 ) 

Extrinsic noise = — = (41b) 

{\ + CV}f 

where the mean mRNA and protein levels are given by (see (39)) 



7 — r k m <^B m y — — k m k p (B,„y , ^-. T a2\ 
<m>= , </>> = (1 + CTJ 

Ym YmYp 
Equation (41a) reveals that when 



(42) 
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<£,„> + ! < 1 

Vm + Vp ~ V 



(43) 



intrinsic noise monotonically decreases with CV^ . On the other 
hand when 



4< J g,„>+l < 1 



(48) 



intrinsic noise first decreases with increasing CF 2 and then 
increases (Figure 4). 



(B m )+l l_ 



Ym + 7/1 "7> 

intrinsic noise first increases with CK 2 , reaches a maximum at 



(44) 



cv:= 



(B„,y + \ 

Ym + Vp 



1 



(B„,) + \ 1 ' 



(45) 



and then decreases with increasing CP" 2 . 

Transcription Burst Size and Translation Rate Fluctuations 

Let transcriptional bursts occur at a constant rate k m with a 
geometrically distributed burst size that is dependent on z(t) and 
given by (21). mRNA translation rate is assumed to be perfectly 
correlated with burst size and is set equal to k p z(t)/(z}. The time 
evolution of moments is obtained from (22) with k p replaced by 
kpz(t)/(z). As in the previous section, although the time 
derivative of fi (Eq. (35)) is not closed, the evolution of fi (Eq. 
(37)) is given by a closed system of linear equations that yield an 
exact expression for the total protein noise level. Recall that 
extrinsic noise is similar for transcription burst size and burst 
frequency fluctuations. Hence, calculation of extrinsic noise for 
model (40) is used to resolve the total noise into its intrinsic and 
extrinsic components. These results show that simultaneous 
transcription burst size/translation rate fluctuations not only 
generate a larger intrinsic noise but also have qualitatively 
different trends compared to burst frequency/translation rate 
fluctuations (Figure 4). 

For slow fluctuations in z(t) compared to mRNA/protein 
turnover rates 



Intrinsic noise = 



<£,„) + ! 

<m> 



1 + 



CV2{l + <B m }(4 + 6CV2 + CV?))\ y p 



(i+(B m y)(i+cvjy 



(46a) 



Ym + Yp 



■ + - 



_ . . CV?(4 + 6CV? + CV?) 

Extrinsic noise = — ~ — , 46b 

(1 + CF?) 2 

where (jri) and (j>y are given by (42). Analysis of (46a) shows that 
when 



4<j? m > + l ; 1 

y m ' yp 

intrinsic noise increases with QVz . However, when 



(47) 



Discussion 

Given the different functional roles of gene expression noise 
inside cells [3,32], much work has focused on understanding how 
variations in the level of a protein arises between otherwise 
identical cells. A class of models were introduced where 
stochasticity arises from two sources: i) Random production and 
degradation of individual mRNA transcripts/protein molecules 
stemming from the inherent probabilistic nature of biochemical 
reactions and ii) Fluctuations in model parameters that correspond 
to randomness in cell-specific factors. Exact analytical formulas for 
total variability in protein level were derived, in spite of the fact 
that in many cases parameter fluctuations lead to nonlinear 
propensity functions. These formulas were decomposed into 
intrinsic and extrinsic noise components as measured by the 
two-color reporter assay (Figure 1). 

Which Mechanism Generates the Largest Gene 
Expression Noise? 

Individual-parameter fluctuations. Comparison of (19), 
(27) and (34) shows that for low values of (,B m }, fluctuations in the 
translation rate create the most variability in protein copy 
numbers. On the other hand for high (-B,,,), burst size fluctuations 
generate the most variability. Burst frequency fluctuations always 
generate the lowest noise. 

Multiple-parameter fluctuations. Equations (41) and (46) 
reveal that simultaneous fluctuations in translation and transcrip- 
tion rates can dramatically increase expression variability. For 
example, consider protein half-life = time-scale of parameter 
fluctuations = 24 hours, mRNA half-life = 8 hours, mean mRNA 
count/cell= 100, <_B m > = 40 and <(p)»<(m). Then, for constant 
parameters, CVj lxed = 0.\ (Eq. (10)). Assuming ATP affects 
transcriptional burst size and translation rate, 10% variability in 
ATP abundance (CV? = 0.1) enhances noise level three-fold from 
0.1 to 0.32. In comparison, burst size fluctuations of similar 
magnitude only increase 0.1 to 0.16. These results reinforce recent 
observations that intercellular variation in ATP abundance can be 
a major driver of gene expression noise [43]. An implicit 
assumption in this analysis is that protein and mRNA degradation 
is insensitive to ATP. Since both ATP-dependent and ATP- 
independent degradation pathways exist within cells, further work 
on ATP-sensitive degradation rates is clearly needed. 

Relationship between Intrinsic Noise and CV^ 

Using Monte Carlo simulation techniques previous studies had 
shown that parameter fluctuations can alter intrinsic noise 
measurements in a two-color assay [42,60]. Building up on these 
results, a systematic analytical analysis of how fluctuations in both 
individual and multiple model parameters affect randomness in 
protein populations counts was performed. Main findings are as 
follows: 

• Intrinsic noise is invariant of fluctuations in the transcription 
burst frequency (i.e., how often mRNA bursts occur from the 
promoter). 
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Intrinsic noise increases with CV Z (extent of parameter 
fluctuations) for fluctuations in the transcription burst size 
(i.e., mean number of mRNAs produced in each burst) or 
mRNA translation rate. 

For simultaneous fluctuations in the burst frequency and 
translation rate, intrinsic noise decreases with CV 1 for low 
protein abundance (Figure 4). Intuitively, for low protein 
abundance (as determined by (43)), the Poissonian term 
1 / mean has a significant contribution to intrinsic noise 
(second term on the right-hand-side of (41a)). Simultaneous 
fluctuations increase mean protein level (see (39)), decreasing 
intrinsic noise. For high protein abundance, ignoring the 
second term in (41a) yields 



Intrinsic noise = 



<£,„> + ! 



1 + 



CV} 



( m ' Yp 



(49) 



which first increases, and then decreases with CV}. The 
maximal value is achieved at CV} = 1. 

Simultaneous fluctuations in the transcription burst size and 
translation rate typically increases intrinsic noise. However, for 
low protein abundance intrinsic noise exhibits a U-shape 
profile with CV 1 (Figure 4). 

In contrast to intrinsic noise, extrinsic noise always monoton- 
ically increases with CV Z . 



CV 



CB m > + l y p 



fixed 



<OT> 

km (But y 

7m 



<p>- 



1 

w 



k p (tri) 



(50) 



Given measurements of CVj ixed and <J?>, a priori knowledge of 
y„ , y m , k p , mean burst size (B m y and frequency k m can be 
computed from (50). Typically, CVj- lxed is assumed to be equal to 
the intrinsic noise measured in a two-color assay. However, our 
results show that this is only valid for transcription burst frequency 
fluctuations. For all other cases, CVj ixed # intrinsic noise, and 
using intrinsic noise for CVj jxed in (50) will lead to erroneous 
parameter estimates [42]. 

Analytical formulas developed here can be used to back 
calculate CVj- lxed from intrinsic and extrinsic noise measurements. 
This point is illustrated for the physiologically relevant parameter 
regime 



<5„,>»1 (Large burst size), (51a) 



</>>»<»?> (High protein abundance), (51b) 



We comment on how these trends change if Fano factor 
(variance/mean), instead of coefficient of variation, is used for 
quantifying noise. This is particularly important in the case of 
multiple-parameter fluctuations, where mean protein levels are 
dependent on CV Z (see (39)). Our analysis shows that in contrast 
to the above trends, the intrinsic noise Fano factor always 
monotonically increases with CV Z for simultaneous fluctuations in 
the transcription and translation rates. 

Recall that our results correspond to a model where mRNAs are 
produced in instantaneous transcriptional bursts. For a promoter 
that stochastically toggles between active and inactive states, this 
approximation corresponds to an unstable active state [47], where 
the promoter quickly transitions back to the inactive state after 
producing a burst of mRNA transcripts form the active state. It 
turns out that some of the above intrinsic noise versus CV] trends 
are also valid outside the instantaneous burst limit. For example, 
Monte Carlo simulations have shown that for fluctuations in the 
translation rate or transcription burst size, intrinsic noise increases 
with CV% when promoter spends a finite amount of time in active 
and inactive states [60]. Future work will extend analytical 
formulas for intrinsic and extrinsic noise to cases where the 
promoter stochastically transitions between different transcription- 
al states. 



CV 2 :« 1 (Small parameter fluctuations). 



In this regime, intrinsic noise is expressed as 



Intrinsic noise = CV% xed {\ +f x Extrinsic noise) 



CVfi xed 'r 



<B m ) Yp 
<m> Yp + 



(51c) 



(52) 



where / = 0 for burst frequency fluctuations and f > 0 in all other 
cases. Analytical expressions for f are provided in the Text SI, and 
it depends only on mRNA, protein turnover rates and time-scale 
of parameter fluctuations (more specifically on ratios Y„/y z and 
Yphm)- Consider a stable reporter protein where y„ = y z = time- 
scale of cell division, and y m = 4y / ,, then 



/=!• 



(53a) 



(Simultaneous fluctuations in burst size and translation rate) 



Estimation of Model Parameters from Noise 
Measurements 

Gene expression noise is often used to calculate the mean 
transcriptional burst size and frequency for a specific gene or 
promoter [40,48,51,52]. Recall from (10) that for fixed model 
parameters 



/=0.4 



(53b) 



(Simultaneous fluctuations in burst frequency and translation rate) 



f =1.6 (Translation rate fluctuations) (53c) 



/ = 2.1 (Transcription burst size fluctuations) (53d) 
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f = 0 (Transcription burst frequency fluctuations). (53e) 



Therefore, if extrinsic noise = 0.5 in an experiment, from (52) 
and (53d), CV^ xed K intrinsic noise/2 for burst size fluctuations. 
Traditional approach of assuming CVj ixed = intrinsic noise would 
overestimate 

CV fixed b y 100% - Using/ = 0.4 for simultaneous 
burst frequency/ translation rate fluctuations gives CVj ixed = 
0.83 x intrinsic noise, and CVj ixed = intrinsic noise may not be 
a bad approximation in this case. It can be shown that 



o</< 



(yp+y m )(Vp+Yz)(yz+y m ) 
y P y m (Vp+Vm+Vz) 



(54) 



with upper (lower) bound being realized for burst size (frequency) 
fluctuations. Without prior knowledge on the source of extrinsic 
noise, (54) yields the following bounds on CVj ixed : 



In conclusion, our analysis reveals how stochastic synthesis and 
degradation of biomolecules combines with parameters fluctua- 
tions to generate heterogeneity in protein level across a clonal cell 
population. These results will help understand how stochastic 
variability is regulated inside cells, and for extracting meaningful 
information from single-cell gene expression measurements. 
Future work will consider scenarios where randomness in cellular 
factor levels simultaneously affects synthesis and degradation 
pathways, or only degradation. Unfortunately, exact solutions are 
unavailable in many of these cases. However, preliminary analysis 
has found moment closure techniques useful for obtaining closed- 
form solutions for the statistical moments. A recent study has 
generalized notions of intrinsic and extrinsic noise from statistical 
moments to temporal correlations [61]. In particular, the auto- 
correlation function of p(i) can be decomposed into intrinsic and 
extrinsic components based on the two-color assay [61]. Future 
will work will derive analytical expressions for protein auto- 
correlation and cross-correlation functions in stochastic models 
with parameter fluctuations, and study how noise signature within 
them can be used for probing genetic systems. 



1 + Extrinsic noise x 



Intrinsic noise 

(v P + y m )(y P + yz)(y z + y m ) 



y P y m {y P +ym+yz) 

< CVp xed < Intrinsic noise, 



(55) 



for the physiologically relevant parameter regime (51). Thus, our 
results provide the necessary correction factors for accurately 
determining CVj ixed from two-color reporter experiments, which 
would be useful for estimating <i? m > and k m . 
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